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Abstract 



We present a complete description of the calculation of the spin-dependent next-to-leading 
order splitting functions. The calculation is performed in the light-cone gauge. We give 
results for different prescriptions for the Dirac matrix 75 in d = 4 — 2e dimensions and 
provide the link to the results in dimensional reduction. 



1 Introduction 



It has become standard to perform analyses of unpolarized parton distributions at next- 
to-leading order (NLO) accuracy of QCD. An indispensable ingredient for such analyses 
are the two-loop splitting functions (or anomalous dimensions) which appear in the NLO 
Q 2 -evolution (GLAP equations. Results (in the MS scheme) for these have originally 

been obtained in ||, using the Operator Product Expansion (OPE) formalism, and in 
where the somewhat more efficient method developed in || was employed^] which is based 
on the factorization properties of mass singularities and on the use of the axial gauge. In 
a recent publication ]§| we have presented a detailed description of the calculation in the 
latter method which had never been fully documented. The main aim of our study was to 
elucidate the role played by the light-cone gauge and the technical aspects related to its use, 
and thus to underline the utility of this gauge for perturbative QCD calculations which had 
been questioned in the past PJlO|]. 

The method of [§|,(|J|] has recently also been applied to derive [|ll| the polarized two-loop 
splitting functions needed for the NLO Q 2 -evolution of the spin-dependent parton densities 
of a longitudinally polarized hadron. Previous OPE results of were confirmed. It is the 
purpose of this paper to provide a more detailed and complete description of our 'polarized' 
calculation jll]] . Although many details of the calculation are the same as for the unpolarized 
case and can therefore be found in Q , there is a new ingredient in the polarized case which 
requires a closer inspection: The Dirac matrix 75 and the totally antisymmetric Levi-Civita 
tensor e^„ p(T enter the calculation as projectors onto definite helicity states of the involved 
longitudinally polarized quarks and gluons. When dimensional regularization is used, a 
prescription for dealing with these (genuinely /owr-dimensional) quantities in d = 4 — 2e 
dimensions has to be adopted which should be free of algebraic inconsistencies. Several such 
'75 schemes', considered to provide a consistent regularization, have been suggested in the 
literature ll^Jl^ , [l5|] . Our calculation [[ll| was performed using the original definitions for 75 
and e^vpa of |l3|l (HVBM scheme) which is usually regarded as the most reliable prescription. 
Besides giving a more detailed account of our previous calculation we will also address the 
use of other 75 schemes such as ]lj,^5| . Furthermore, we will provide the link to the results 
in dimensional reduction. 



2 The Calculation 

2.1 Framework 

An outline of the method of |],|||| to calculate (NLO) splitting functions as well as a detailed 
description of the calculation in the unpolarized case have recently been given in ||. We 
will therefore only focus on the new aspects arising in the polarized case, some of which 
have already been discussed in our previous paper . 

As usual, all polarized quantities like cross sections etc. will be denoted by a 'A', i.e., 

AM = i (M(++) + M(— ) - M (+-) - M (-+)) , (1) 

1 Note that thereis a discrepancy between the results of |J and which was resolved in jj] in favor of 
the calculation of [H. 
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where '+, — ' stand for the helicities of the scattering incoming particles. The polarized 
parton distributions A/ (/ = q, g) are defined by 

A/ = /+-/- , (2) 

/ + (/ — ) denoting the density of parton-type / with positive (negative) helicity in a nucleon 
with positive helicity. Omitting the 'A' in Eqs. (|l]j2|) and taking the sum on the right-hand- 
sides, one recovers the analogous relations for the unpolarized cross sections and parton 
distributions. 

The general strategy consists of first expanding the squared matrix element AM for 
polarized virtual photon-polarized quark (gluon) scattering into a ladder of two-particle 
irreducible (2PI) kernels B Co, K , 



AM = A 



C a {l + K G + Kl + Kl + ...) 



EE A 


Co 


1-Kq 





We now choose the light-cone gauge by introducing a light-like vector n (n 2 
n ■ A = 0. At the same time, n is used to define the longitudinal direction: 

n ■ p ee pn 0, n ■ t = p ■ t = , 



(3) 



= 0) with 



(4) 



where p is the momentum of the incoming parton (taken to be massless), and t is any 
vector in the transverse plane. In the light-cone gauge the 2PI kernels are finite before the 
integration over the sides of the ladder is performed. Collinear singularities therefore appear 
only when integrating over the lines connecting the rungs of the ladder || . This allows for 
projecting out the singularities by introducing the projector onto polarized physical states, 
AT. More precisely, AT decouples the product, A(AB), of two successive 2PI kernels by 
projecting onto definite helicity states of the particle connecting the kernels and by setting 
this particle on-shell in A. Writing down explicitly the combinations of the helicities of the 
in- and outgoing partons and of the intermediate particle, one immediately obtains from 
Eq. ([!]) A(AAP B) — > AAAB, i.e., the (decoupled) product of two polarized kernels. Thus 
AM can be written in the factorized form 



AM = ACAr , 

where (introducing the modified kernel K — Kq(1 — (1 — AP)Kq)^ 1 ) 

1 



AC 

Ar 



ACo 



1 - (1 - AV)K 



1 



1 - AVK 

1 + AVK + AVK {1 - AV)K + (APK )(APK ) 



(5) 



(6) 



(7) 



AC is interpreted as the (finite) short-distance cross section, whereas Ar contains all (and 
only) mass singularities. Working in dimensional regularization (d = 4 — 2e) in the MS 
scheme one has explicitly: 



ATij(x,a s , -) = Zj 



(5(1 - x)8 lJ +jPP 



d d k 
(2n)° 



5{x 



n ■ k , 
pn 



)AU l K 



l-AVK ALl 



, (8) 
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where 'PP' extracts the pole part of the expression on its right and Zj (j — q(g)) is the 
residue of the pole of the full quark (gluon) propagator, k is the momentum of the par- 
ton leaving the uppermost kernel in Ar; by definition of n, x can be interpreted as the 
infinite-momentum frame (IMF) momentum fraction of p carried by k. The spin-dependent 
projection operators onto physical states are given by 

AC/g = ~^~7^75^, AL 9 = -2% 

AU g = ie^^, AL g =ie^^ . (9) 
n-k 2pn 

We see that the quantities 75 and e Atl/pCT appear, for which we will have to define a continu- 
ation to d dimensions. 

Finally, it can be shown Q that the coefficient of the 1/e pole of Ar corresponds to the 
GLAP [|]J|] evolution kernels we are looking for: 

AT qq {x, a s ,e) = 5(1 - x) -\ (^-AP^(x) + \ (g)' AP$(x) + . . + O , (10) 

and analogously for the flavor singlet case. Here we have adopted the perturbative expansion 
of the splitting functions, 

AP l3 {x, as )=(^)AP^\x) + (^-) 2 AP^(x) + ... . (11) 

We conclude this section by collecting all ingredients for a NLO study of longitudi- 
nally polarized deep-inelastic scattering in terms of the spin-dependent structure function 
gi(x,Q 2 ). There are two different short-distance cross sections^, AC q and AC g , for scat- 
tering off incoming polarized quarks and gluons, respectively. They are to be calculated 
according to Eq. (0). Thus gi reads to next-to-leading order: 



+ 



1 n/ 

-Y 

i=l 

a s (Q 2 ) 



9l (x,Q 2 ) = ty j( % { Ag^Q^ + A^Q 2 ) 



i=l 

2\ r 



AC„ <g> (Aq, + Aq t ) + —AC q ® Ag 



(x,Q 2 )\, (12) 



2- 

where iif is the number of flavors and <g> denotes the convolution 

f<Z>g= [ dydz f(y)g(z) S(x - yz) . (13) 
Jo 

Defining the sum and the difference of polarized quark and antiquark distributions as 

Aqf = Aq % ± Aqi , (14) 

one finds the following evolution equations for the non-singlets Aq~ and Aqf — Aq^ (see, 
e.g., 1): 

d 



d\nQ< 



(Aq+-Aq+)(x,Q 2 ) = AP q + q (x ias (Q 2 ))®(Aq+ ~Aq+)(x,Q 2 ) , (15) 
'' Aqr{x,Q 2 ) = AP q -(x.a s {Q 2 ))® Aqr(x,Q 2 ) , (16) 



dlnQ 



"We do not distinguish between a quark non-singlet and a quark singlet short-distance cross section here. 
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where 

AP± = APV ± AP£ , (17) 

with AP^ q starting to be different from zero beyond the leading order. Introducing the 
polarized quark singlet AE = J2i(^Qi + Aq~i) one nas m the singlet sector: 




The qq entry in the singlet matrix of splitting functions is given by[] 

AP qq = AP+ + AP q s q . (19) 

So, at NLO, we will have to derive the splitting functions AP^'^, AP^ 1 ^, and those 
involving gluons. 



2.2 NLO graphs 

According to Eqs. (0J|), Ar is given to NLO by 

Ar = Zj (1 + APA" + AP(A 2 ) - AP(A APA )) . (20) 

The basic topologies of all 2PI diagrams which occur in NLO are shown in Fig. 0, where 
the notation of the topologies is as in J8| . Explicit examples of graphs contributing to the 
various splitting functions can be found in @Jll]|| . Topologies (hi) correspond to the terms 
AP(A^) - AP(A APA ) in Eq. @, all other topologies belong to APA . As is obvious, 
topologies (cd) and (fg) possess a real and a virtual cut. Fig. [l] does not display the genuine 
two-loop graphs which determine Zj and thus the endpoint (5(1 — x)) contributions to the 
diagonal splitting functions. 




(b) (cd) (e) (fg) 




(hi) (jk) 
Figure 1 : Basic topologies of the NLO diagrams 



3 As compared to |j we include a factor 2nj in the definition of AP C 
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2.3 Use of the light-cone gauge 

As noted in the previous section, the light-cone gauge plays a crucial role in the calculation. 
In this gauge the gluon propagator takes the form 

where I is the gluon's momentum. As is well-known, the light-cone gauge propagator 
can give rise to additional divergencies in loop and phase space integrals. We follow 
to use the principal value (PV) prescription to regulate such poles: 

1 1 \ _ n ■ I 

n ■ I + i5(pn) n ■ I — i6(pri) J (n ■ I) 2 + S 2 (pn) 2 

All singularities related to the gauge propagator can then be cast into the basic integrals 

§1 

= (-0,1). (23) 

The PV prescription leads to the feature that the renormalization constants depend on Iq 
and the longitudinal momentum fractions x The complete set of the renormalization 

constants, as well as most other technical ingredients we need for our calculation, like phase 
space integrals and scalar virtual integrals with and without a gauge propagator, have been 
worked out and collected in 

2.4 Treatment of 75 and e pupa 

As seen from Eq. (^, the crucial difference with respect to the unpolarized calculation is 
the presence of 75 and t^upa which act as projectors. Definitions for these quantities in 
d = 4 — 2e dimensions (or at least rules for handling them) have to be given. 

It is well known that a 'naive' fully anticommuting (i-dimensional 75 can easily lead to 
algebraic inconsistencies when located in a trace with six or more other Dirac matrices | fl6| , 
and thus results obtained this way cannot really be trusted. The problem can be avoided in 
essentially two ways: 

• One can define 75 in d dimensions by maintaining its /owr-dimensional definition 

75 ee le^ JlM j u j p7a , (24) 

where the e-tensor is regarded as a genuinely four-dimensional object, i.e., its com- 
ponents vanish in all unphysical dimensions. This definition of 75 is the original one 
of [|l3| (HVBM scheme). Splitting the d-dimensional metric tensor into its four- and 
(d — 4)-dimensional components, 

9ij.u = 9 pv + 9^ (where = 4, gfi = d - 4) , (25) 

and defining 7^ = g pv ~i v etc., one finds from ( pi| ) that 

{7^,75} = for n = 0,1, 2, 3, 
[7", 75] = otherwise. (26) 
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• One can maintain the full anticommutativity of 75 with all other Dirac matrices but 
abandon cyclicity of the trace operation iQ. In this case a fixed 'reading point' has 
to be defined from which all Dirac traces corresponding to a given process have to 
be started. As in the HVBM scheme, the e-tensor is defined as an entirely four- 
dimensional object p"4j . 

Our main calculation 4 is done in the HVBM scheme since its consistency is well-established. 
We will also discuss the results obtained in the scheme of JlJ] which was used in the OPE 
calculation of the AP^ 1 ' of fL^| . Due to the fact that the Levi-Civita tensor is assumed 
to be four-dimensional (and because of ( |26| ) in the HVBM scheme), the matrix element 
of a graph will for both prescriptions not only depend on the usual 'd-dimensional' scalar 
products of two momenta, like I ■ k = g^vl^k 11 etc., but also on '{d — 4)-dimensional' ones, 
denoted by I ■ k = g^l^W ', k etc. Special care has to be taken to take into account 
such (d — 4)-dimensional terms in loop and phase space integrals. Details are given in the 
appendices. 

We mention that in Jll],[l1| another prescription was suggested to handle traces with one 
75. One eliminates 75 via the relation 

7^75 = 7^ W7V7" • (27) 

The remaining trace is then perfectly well-defined. One always ends up with the product 
of two Levi-Civita tensors which can be written as a determinant of a matrix, the elements 
of which are metric tensors. Following |jl5||lS| ], one then performs the contractions of these 
metric tensors in d dimensions^. The result should then coincide with that in the HVBM 
scheme. The scheme is attractive because it avoids any (d — 4)-dimensional scalar products, 
but it leads to longer trace operations because of the substitution (p7|). We will briefly 
return to this prescription when presenting our results. 

We finally note that a variant of dimensional regularization, dimensional reduction [ fl9[ , 
has been widely discussed and used in the past years, the main reason being its applica- 
bility to supersymmetry. The scheme essentially consists of performing the Dirac-algebra 
in four dimensions, which makes the treatment of 75 straightforward, and of continuing 
only momenta to d {d < 4) dimensions. In order to match the ultraviolet (UV) sectors 
of dimensional regularization and dimensional reduction, specific counterterms have to be 
introduced |^0|,|l| in the latter. Once this is done, there is a straightforward and universal 
way to deal with differences arising from mass singularities [|l 22 23j. We will exploit the 



results of ]21|,^2j,g3[ in order to translate our results to the form they take in dimensional 
reduction. 



3 Results 

The full one-loop results are included for completeness ]l]j24| : 

AP^(x) = C F {^-r--l-x+~6(l-x)} (28) 

4 We use the program Tracer of ||l7{ for calculating the Dirac traces and .performing contractions. 
5 Here one has to multiply the gluonic projection operator AL g in feh by a normalization factor 
2/((d-2)(<f-3)). 



G 



where 



a^ 

A^ 0) ( a 

AP(°)(a 



27> {2x - 1} 

C F {2 - x) 
1 



= 2V r 



2x + l} + tj5(l - x) 



Cf = n c = 3, T> = T R n f = |n/, /3 
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(29) 
(30) 

(31) 
(32) 



3.1 Results in the HVBM scheme 

The graph-by-graph results (for those splitting functions that involve more than just one 
topology) in the HVBM scheme are given in columns (b)-(jk) of Tables 1-4. For the non- 
singlet case, Table 1, we only present the differences between our graph- by-graph results 



for the polarized APY q ^ and the corresponding contributions to the unpolarized Pqq^ 1 as 



listed in Table 1 of g . The reason for this is that for an anticommuting 75 the qq diagrams 
would trivially yield AP^ 1 ' 



Pqq^ because of the fact that there are always two 75 



in the same trace which could be anticommuted towards each other and eliminated using 
7g = 1. However, as discussed in section 2.4, 75 does not fully anticommute in the HVBM 
scheme. So it is interesting to see what happens to the difference AP L 
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P 



v,(i) 
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when this 

prescription is used. To examine this question further we have also distinguished in Table 
1 the contributions coming from the [d — 4)-dimensional scalar products, like k etc. (see 
section 2.4 and the appendices), from all other contributions^. 

Let us concentrate on the columns 'Sum' in each table, which give the sums of the 
graph- by-graph results (b)-(jk). The first thing to observe is that all contributions from 
the integrals Iq, 1%, arising from the light-cone gauge propagator, cancel in 'Sum', as they 
must. We furthermore see from Table 1 that APj£^ — Pqq is indeed non-vanishing in 
the terms multiplied by CfNc and CfT/. To study the implications of this let us first 
note that the endpoint contributions (~ 8(1 — x)) to the diagonal splitting functions are 
necessarily the same as in the unpolarized case (where they were derived in [p5p|]), since 
they are determined by Zj in Eq. (^). We can then express the results corresponding to the 
columns 'Sum' in terms of the unpolarized NLO non-singlet splitting functions P^q of 0] 



and the recent polarized OPE results ApV of 



AP±/V(x) 


= PT q {1 \x)~2foC F (l-x) , 


AP q s q M(x) 


= APf/ 1 ^) , 


AP«(z) 


= AP« (a;) + 4C F (1 - x) ® APW (x 


AP«(z) 


= AP« (x) - 4C F (1 - x) ® APW (x 


APW(x) 


= AP«(x) , 


ACq(x) 


= ACq(x) - 4C F (1 - x) , 


AC g (x) 


= AC 9 (z) , 



6 In all other tables the contribution s from the (d — 4)-dimensional scalar products are aire 
7 In our normalization given by Eq. (hllha) the results for the NLO splitting functions of t 
divided by 8 and those for the NLO short-distance cross sections AC q , AC 9 by a s /2n. 



(33) 
(34) 
(35) 
(36) 
(37) 
(38) 
(39) 

iy included. 
H have to be 



7 



where we have also included the results for the short-distance cross sections AC q , AC g 
which are to be calculated according to Eq. (||). Obviously, the term — 2/3oCf(1 — x) in 
( |33| ) is entirely due to the fact the HVBM 75 does not fully anticommute. As was already 
discussed in p7| , p8[ and indicated in Eq. (p3|), the '+' and '— ' combinations of the NS 
splitting functions as defined in ([17]) interchange their role when going from the unpolarized 
to the polarized case, equivalent to AP^'^ = — Pj* . The latter relation is a consequence 
of the projection operator for antiquaries j|, +75 (— fCj/in ■ k, in this case. Again it is trivial 
for a fully anticommuting 75, but it turns out that it is also respected by the HVBM 75. 
Eqs. therefore imply that the combination AP^ {1) = P w ' (1) - 2/3 C F (l - x) 

would govern the Q 2 -evolution of, e.g., the polarized NS quark combination 

AA 3 (x,Q 2 ) = (Au+ - Ad+) {x,Q 2 ) . 

Since the first moment (i.e., the s-integral) of the latter corresponds to the nucleon ma- 
trix element of the NS axial vector current q^y^^^X^q which is conserved, it has to be Q 2 - 
independent [ffijf . Keeping in mind that the integral of the unpolarized Pqq'^ vanishes 
already due to fermion number conservation Q, it becomes obvious that the additional 
term — 2/3oCf(1 — x) in (|33| ) spoils the Q 2 -independence of the first moment of AA 3 (x, Q 2 )- 
It is therefore necessary to perform a factorization scheme transformation to the results in 



(p3-39) in order to remove this additional term. Such a scheme transformation is always 



allowed since neither the AP^ nor the ACi are physical quantities. Thus one can shift 
terms between them in a well-defined way without changing a physical quantity like gx, 
hereby just redefining the polarized NLO parton distributions. Even though it is a priori 
only the term — 2/3 Cf(1 — x) in ( |33| ) we want to remove, the transformation will also affect 
the singlet sector since, according to Eq. (|l9|), APqq 1 ^ also occurs in the singlet evolution 
matrix. The scheme transformation reads in general (see, e.g., [p0||l2|]): 



qq 


AP±'« - 2p Az qq 


1 










APW - 
qq 


AP$-2p Azqq + 


AAZqg d 


, APW - 


4Az 3? (g 


AP(°) , 




APW - 
qg 


APW - 2P AZqq + 


AAZqg d 


5 (APW 


- APi?; 






AP« - 

gg 


- AP£ -2f3 Q Az gq + 


AAZgq (> 


s (ap%> 


-ap(°); 




5 (Az ff9 


AP« - 

gg 


- AP« - 2(3 Az gg - 


AAZqg i 


5 AP^ + 


4Az gq (8 


Ap(°) 




AC q - 


-► ACq ~ AAZqq , 












AC g 


-» ACg - iAZqg , 













AZgg) , 

AZgg) , 

(40) 

where the Azy generate the transformation. For the case at hand we only need a non- 
vanishing AZqq, 

AZqq(x) = -C F (l - X) . (41) 

The changes to our result caused by inserting Az qq into (^) are given in the columns "75" 
in Tables 1-3, to be added to columns 'Sum' to obtain the final answer. As one can see from 
Table 1, the difference APjg ' —Pqq finally becomes zero, which is a result of cancellations 
between terms from (d — 4)-dimcnsional scalar products and other terms. Furthermore, it 
turns out that the transformation ( |40| , [4l| ) not only removes the term — 2/?oCf(1 — x) from 
Eq. (|33|), but eliminates all extra (1 — x)-terms on the r.h.s. of (p3[]39|), leaving AP^'*- 1 ^, 
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APpg^ and AC g unchanged. Thus our final results are in complete agreement with those of 
|p^| . We finally note that the presence of the (1 — x)-terms in our original HVBM scheme 
result (R3m9) can be traced back to the fact that in this scheme the polarized LO splitting 



(0),<J=4-2< 



function in d — 4 — 2e dimensions, AP qq 
counterpart, i.e., violates helicity conservation: 



, is no longer equal to its unpolarized 



AP (0)4^-2e {x) _ p(0),^4- 2%T) = 4e(1 _ x) 



(42) 



The additional term — 4(7f(1 — x) in the HVBM-scheme result for AC q in ( |38| ) was already 
identified in pl],p8[ . After its removal by the transformation ([h^,|4]), the integral over AC q 
takes the value — 3Cf/2, giving rise to the correct NLO correction (1 — a s /n) to, e.g., the 
Bj0rken sum rule. 

Our complete final results can now be collected from columns 'final' (or 'Sum' if there 
has been no change due to (|I],|d|)) in Tables 1-4. To write them down we introduce 



Sp q g(x) = 

Spgg(x) = 

Sp gg {x) = 



2x-l 
2-x 
1 



- 2x + 1 



(43) 



We then have 



AP±.(1) = PT,(1) 
99 99 



APW\x) = 2C F T f 



APi^x) 



(1-a;) - (1- 3x) lux - (l + x) In 2 x 



22 + 27x - 9 In a; + 8 (1 - x) ln(l - x) 



C F T f 



+Sp qg (x) 2 In (1 - x) - 41n(l - x) In x + In x tt' 



+N c T f 



2 (12 - llx) - 8 (1 - x) ln(l - x) + 2 (1 + 8x) \nx 



-2 ^ln 2 (l - x) T —j Sp qg (x) - {2S 2 {x) - 3 In 2 x) Sp qg (~x) 



C F T f 



+Cl 



-(x + i)--6p gq {x) ln(l-x) 



1 1 



(4 — x) In x — 5p gq (—x) ln(l — x) 



2 2 

j i-\vL 2 {\ - x) + ^\n 2 x) 5p gq (x) 



+C F N C 



1 



(4 - 13a;) lnx + i (10 + x) hi(l - x) + ^ (41 + 35x) 



+ -(-2S 2 (x)+3\n 2 x)5 Pgq (- 



(44) 
(45) 



(46) 
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In (1 — x) — 21n(l — x) In a: — j Sp gq (x) 



(47) 



AP$(x) = -N c Tf 



-C F T, 



4 20 4 

4(1 - x) + - (1 + x)lnx + -y5 Pgg {x) + -6(1 - x) 



10 (1 - x) + 2 (5 - x) Iilx + 2 (1 + x) In 2 x + 6(1 - x) 



c 
G7 



\ (29 - 67a;) In a; - ^ (1 - x) + 4 (1 + x) In 2 x - 2S 2 (x)5p gg (-x) 



4 ln(l — x) In x 4- In x 

9 3 



djp OT (aO+ 3C(3) + - )6(l-x) 



AC q (x) = C f 



(1 + x 2 ) 
+ 2 + x 



ln(l - cc) 



1- a; 



1 + x 2 



2 [1 - x]+ 1 - x 



In a; 



AC s (a;) 



27V 



(2a; - 1) ( In - — - - 1 ] + 2(1 - x] 



where the unpolarized NS pieces Pq\ can be found in 

S 2 (x) 

For relating the above results to those of |12j the relation 



and 



1 +* dz .1 — Zs 
— ln( ) 

a Z Z ' 

1+X 



1 , 2 rr- 



S 2 (x) = -2Li 2 (-a;) - 21na;ln(l + x) + - In x ■ ( . 



(48) 

(49) 
(50) 

(51) 
(52) 



is needed, where Li2(a;) is the Dilogarithm P6[ . As explained above, the 5(1 — a;)-endpoint 
contributions could be taken from the unpolarized case f|[]||; £(3) ~ 1.202057. The +- 
prescription in ([13049) is defined in the usual way, 



dzf(z)[g(z)},= / dz(f(z)-f(l))g(z) 



(53) 



It is obviously only needed if the function multiplying it is non- vanishing at x = 1. 

For completeness we finally list the first moments 8 (x-integrals) of the results in Eqs. (E5| 

J AP^(x)dx = -3C F 7> , J AP£>(x)dx = ~lc 2 F + 7 -±N c C F - l -C F T f , 



8 Compact expressions for all Mellin-moments of the NLO splitting functions, J dxx 11 'AP^'fi), and 
their analytic continuations, can be found in p7[ . 
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J 1 APg\x)dx = , J 1 AP$(x)dx = l ^Nl ~ C F T f - b -N c T f = ^ , 

f AC q (x)dx = ~C F , f AC g {x)dx = 0. (54) 
Jo 2 Jo 

3.2 Results for other 75 prescriptions 

We now discuss the results we obtain when using the prescription of ]l4| with an anticom- 
muting 75 and a non-cyclic trace. As mentioned earlier, the property {7^,75} = will 
automatically yield AP^'^ — P^j instead of (|33|). Thus it immediately follows that 
either all other (1 — a;)-terms in ( ^3[]39| ) are absent as well, or that the full result will be 
genuinely different from the one in Eqs. (|44] - ^o|) , i.e., not transformable into (|4^ - ^o|) by a 



factorization scheme transformation (40). To study this question, we have to calculate the 



other NLO splitting functions and the short-distance cross sections. Since the prescription 
for the Levi-Civita tensor is the same as in the HVBM scheme, it is obvious that APgg 
and AC g will be the same in both schemes. The remaining quantities are the interesting 
ones since they involve one or two traces with one 75. In the prescription of |l4[ , one has 
to define a 'reading point' for these, at which the trace is to be started. Let us discuss our 
choice for the case of the C F Tf part of AP q l\ in which the incoming particle is a gluon and 
the outgoing one a quark with the projector —75^/471 • k being acted upon. Only topologies 
(cd),(e),(fg) and (hi) contribute; for the sake of clarity we show the graphs explicitly in 
Fig. |^. It seems reasonable to choose the projector as the reading point, marked by the 




'crosses' ffi in Fig. ^. Doing so for graphs (cd),(e),(fg), which in the language of Eq. ( |20| ) 
represent AVKq, one obtains for each graph exactly the same answer as when using the 
HVBM 75. However, for topology (hi) one has to take care: For the subtraction graph (i), 
which represents — A'P(KqA'PKq), the projection operator AV has by definition effectively 
cut the Dirac-trace of (h) (= A'P(Kq)) into two traces, hereby inserting an additional 75 
in each trace according to AL q , AU q in ([)[). Thus (i) essentially becomes the convolution 

of APg°' for the upper part of the diagram with APg!, for the lower. Since the trace for 
the upper part then contains two 75 which eliminate each other, the 75-problem has been 
shifted from the upper 'cross' to the lower one. This suggests that the trace for graph (h) 
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should be read from the black dot rather than from the 'cross'. Indeed, when doing this, 
it turns out that both pieces (h),(i) yield differences with respect to the HVBM results, 
which on aggregate are exactly the same as the entries in the column "75" in Table 2. This 
means that we end up with the result in dlrjj ) for the CpNf part of APqg , without any 
extra terms. The same thing happens for the case of the C\ part of APgq , where it is 
resonable to choose the —3^75 of the incoming polarized quark as the reading point for the 
graphs for AVKq, in which case each of these graphs individually gives the same result as 
in the HVBM scheme. This time, the answer for subtraction graph (i) is also the same as 
in the HVBM scheme, but when the reading point for graph (h) (AV(Kq)) is again chosen 
as the quark propagator on the side of the ladder between the rungs, it alone generates all 
differences in column "75" of Table 3, such that the C\ term in ( |47| ) is reproduced. Next, 
we have to calculate the NcTf part of AP^ and the CfNc part of APgl\ choosing again 
— 75^/4n-/c and —^75 as the reading points, respectively. Here, no extra care has to be taken 
since for both the sides of the ladder between the rungs in graph (h) are gluons. For each 

graph the result exactly agrees with the corresponding one in the HVBM scheme. Finally, 

s (1) 

we have calculated AP qq , again choosing the projectors as the reading points for its two 
traces and reproducing once again the result of the calculation in the HVBM scheme. Also, 
AC q in ( [49| ) was obtained in [^8| with the photon vertex as reading point. 

Thus, to summarize our results for the calculation in the 75-scheme of |l4j], we have 
indeed reproduced all polarized NLO quantities as given in their final form in Eqs. (fi~i|-f5o|) . 
without any extra terms like in d33|-[39|) . This is a nice confirmation of our HVBM results 
and underlines the consistency of the whole calculation. However, it needs to be emphasized 
that although our choice of the reading points above seems the most reasonable one, other 
choices should be allowed. But if one insists on choosing the reading point — 75^/471 • k also 
for graph (h) of the CpTj part of AP^ (rather than shifting it to the black dot in Fig. ||), 
the result for this graph is no longer different from the one in the HVBM scheme! Thus a 
genuinely different final answer for APq^ would be the consequence. We also found that 
choosing a different reading point in the calculation of AC q the result seems to change. It 
appears likely that these problems indicate a certain incompatibility of the 75-scheme of |lj] 
with the projection method we are usingP] in which, as discussed above, traces sometimes 
are cut, with insertions of 75, by action of AV . The genuinely four-dimensional (i.e., non- 
anticommuting) 75 of |l5| ] seems to be more appropriate here since it gives a unique answer 
without any extra effort apart from the scheme transformation (|4^^l|) . 

We mention that we have also recalculated some graphs of the CfTj part of APqg , as 
well as the coefficient function AC q , using the prescription of flli^Jlq], i.e., eliminating 75 via 



Eq. (27) and contracting the Levi-Civita tensors in d dimensions. This scheme is expected to 
be equivalent to the HVBM scheme |l5|] . In all cases we studied we indeed obtained the same 
answer as in the HVBM scheme without having to take into account any (d — 4)-dimensional 
scalar products which are not present in this scheme since the metric tensors coming from 
the products of two e-tensors are taken to be d-dimensional. This computational advantage 
is lost, however, because the trace operations become considerably more involved due to the 
substitution (|27|). Furthermore, it is not entirely clear how to generalize this scheme in order 
to deal with traces with two 75, i.e., how to establish equivalence to the HVBM prescription 
in this case. 



9 On the other hand, similar problems with the prescription of [|L4j also seem to have been the source of 
an error in a previous version of the OPE calculation of pa]. 
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3.3 Connection with dimensional reduction 



We finally discuss how our results in Eqs. fl44| - |5C| ) can be translated to dimensional reduction. 
As was explained in 20,2]]], the UV sectors of QCD in dimensional regularization and 
dimensional reduction are made to agree by introducing additional counterterms in the 
latter which include a finite renormalization of the strong charge. Once this is done, all 
remaining differences between the results for a NLO quantity in dimensional regularization 
and in dimensional reduction can only be due to the effects of mass singularities. They 
are fully accounted for ^l],^2|,^3| by the differences between the d-dimensional LO splitting 
functions (as to be obtained in dimensional regularization) and the /our-dimensional ones 
(corresponding to dimensional reduction). This makes it very easy to transform our results 
in Eqs. ( p^ - |50| ) to dimensional reduction: We just need to perform a factorization scheme 
transformation ([h]), with the Az%j to be obtained from the parts ~ e of the polarized 
d-dimensional LO splitting functions as obtained in the HVBM scheme. The latter read: 

A p(0U^-2e {x) = CF {_Jl__ 1 _ x + 3e{1 _ x)+ ^l S{1 _ x) } (55) 

AP (0),d=4-2e {x) = 27>{2x- 1-26(1-*)} (56) 
AP(;)' d=4 - 2( (i) = C F {2 - x + 2e(l - x)} (57) 

A p(0),<*=4-2e W = 2Nc ^_l__ 2x + 1 + 2e{1 ^ x) ^ + ^ + ^ £ y {1 _ x)m 



However, care has to be taken: We have already performed the scheme transformation ( |4l| ) 



which effectively amounted to removing 4e(l — x) from APq^' d 4 2c , hereby rendering it 



equal to its unpolarized counterpart, Pq q r> ' d 4 2e (see Eq. (|42|)). We thus have to use 

AP (0), rf =4-2c (a;) = C i ,{ ir ^- r -l-x-e(l-x) + *±±8(1 x)} (59) 



instead of (|55|). From the parts ~ e in Eqs. (|56|-p9D we can now read off the Az^ to be used 
in the new scheme transformation (|40|): 

Az m (x) = ^((l-x)-±S(l-x)\ , Az qg (x) =T f (l-x) , 



Az gq (x) = -^(l-x) , Az gg (x) = -N c ((l-x) + ±6(l-,)) . (i.m 



Inserting these into (pl)[), one finds: 

A p(l),DR„ Ap (l) (x) = _^^ 1 _ x) _l S{1 _ x ^_ 4CFTf{1 _ x)hxXj (61) 

A pW,DR_ &pW( x ) = C F Tf (9 - Wx + 6 lna;- 8(1- x) ln(l - x)) 

/ 25 26 \ 
+ N c T f (-— + —x + 8(1- x)(hi(l-x)- lnx)\ , (62) 

AP^ q ) BR - AP S (1) (x) = C F U - -x - 4(1 - x) ln(l - x) + (4 - x) In x 
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APi 



1),DR 



.•/;/ 



/ 19 \ 8 

+ C F N C f 3 - —x + 4(1 - x) ln(l - x)\ - ~C F T f (l - x) (63) 

ApW(.t) - 2foN c ({l-x) + ^8(l-x)\+4C F Tf(l-x)ln.x, (64) 



where the AP^- are our final results in Eqs. (|44|-[48|) and the AP^p' DR are the polarized 
NLO (MS) splitting functions in dimensional reduction. For the qq sector we have only 
written down the entry of the singlet matrix, AP$ = AP+' (1) + AP<f g ' (1) (see Eq. @). 
Furthermore, we have not written down the short-distance cross sections since their transfor- 
mation is trivial. We mention however, that the effect of Az qg on AC g is to remove the term 



+2(1 — x) at the end of (fjOj), which corresponds to a factorization scheme in which gluons 
contribute to the first moment (x- integral) of the structure function g\ . Such a factorization 
scheme was suggested in [^2|. Since furthermore the first moments of APqq- > ' DR and APq^' DR 
vanish, it follows that the total polarization of quarks and antiquarks, dxAT, DR (x, Q 2 ), 
is Q 2 -independent when defined in dimensional reduction. 

As a little cross-check on our findings in Eqs. (^3l||64|) we have recalculated all polarized 
NLO splitting functions, using all matrix elements in four dimensions (i.e., setting e and (d— 
4)-dimensional scalar products to zero), but performing loop and phase space integrations 
in d dimensions after all spin algebra has been done. Subtracting the corresponding results 
from our HVBM ones, Eqs. (Q-^8|), we reproduced all logarithmic terms in Eqs. (|6l[-p4|), 
leaving only some non-logarithmic pieces (and of course the ^-functions) unaccounted for. 
It is precisely such non-logarithmic terms one expects to be generated by inclusion of the 
additional finite UV-counterterms needed for a full and proper calculation in dimensional 
reduction MM. 



The main reason for performing our exercise concerning dimensional reduction is that 



the results for AP^' DR to be obtained from Eqs. (61 - 64|) , (|44-|48|) fulfil a remarkable relation 



if one sets C F = N c = 2T f = N (see also [L2 



AP$>™(x) + AP^(x) - AP^(x) - AP^(x) (65) 

which is known to hold for LO splitting functions and also for the unpolarized NLO splitting 
functions in dimensional reduction J33|,[22[ , and is expected from supersymmctry. 



4 Summary 

We have presented a detailed description of our calculation [|ll| of the spin-dependent NLO 
GLAP splitting functions in the method of |||]|^] • The main calculation has been performed 
using the 75-prescription of p3| , but we have also discussed the results one obtains when 
using a fully anticommuting 75, giving up cyclicity of the Dirac-trace |jl4|| . The results 
for the two schemes turn out to be the same for a reasonably chosen reading point in the 
latter scheme, but the HVBM scheme appears to be the safer and the more straightforward 
prescription. Our final result confirms the previous OPE results of Jl^| . We have also 
provided the connection with dimensional reduction, in which the NLO splitting functions 
satisfy a simple supersymmetric relation. 
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Appendix 



Almost all technical ingredients we need, like two- and three-point functions and three- 
particle phase space integrals, have recently been given in the documentation M of the 
unpolarized calculation. We only need to consider some new details related to the treatment 
of 75 and e'" yp<T in d dimensions, which can introduce explicit dependence of the matrix 
element squared on scalar products of the nonphysical (d — 4-dimcnsional) components of 
vectors (see section 2.4). It is most convenient to work in the IMF parametrization of the 
momenta 01: 



P 
n 



( P , xy , 
pn 




-, k 



where p represents the incoming and k the outgoing ('observed') parton, with 

k 2 = k 2 x + k 2 y + {kf = k 2 T + k 2 



(A.l) 



(A.2) 



being the total transverse momentum squared of k relative to the longitudinal axis defined 
by p, n. Here we have explicitly introduced the (d — 4)-dimensional components of k with 



{k) 2 = -g^k v , 



(A.3) 



9tiu being the (d — 4)-dimensional metric tensor of section 2.4. The labelling of the other 
momenta is fixed in Fig. [|. 




p-k 




Pa 




(a) (b) (c) 

Figure 3: (a) Vertex correction graph (b) One parton emission (c) Two parton emission. 



A Virtual Integrals 

The virtual diagrams always comprise emission of only one real massless particle (the rung 
of the ladder) from the process p — > k + p^ (see Fig. ||(a),(b)). The corresponding phase 
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space reads (including the integration over the 'observed' parton k): 

PSW = J d d k x5(x -^)J j^h5 + {pl) {^) d 5 d { P P s k) 

' (l-x)- e [ Q d\k 2 \\k 2 \-4(-e) [ dkk- 1 -^ , (A.4) 



r(i-e) 

where h is defined by 

k 2 = \k 2 \(l -x)h . (A.5) 

The last integral in Eq. (A.4) has been written in such a way that it is unity if there is 
no dependence on k, i.e., k . As we can also see, any term proportional k in the matrix 
element will give a factor e after performing the K-integration. 

When performing a loop calculation in the polarized case, the matrix element will in 
general not only depend on k 2 , but also on scalar products involving the (d — 4)-dimensional 

components of the loop momentum i, like I 2 , k ■ I. One can always write, e.g., I ■ k as 
—g finish" , then perform the loop integration, and contract the result with —g l _ iv k v after- 
wards. Obviously, this requires knowledge of vectorial two- or three-point functions (and 
tensorial ones when I 2 = —g iiv l li l 1 ' appears). These can be easily obtained from the scalar 
ones presented in The only new feature is that, unlike in the unpolarized calculation, 
we also need the vectorial three-point functions with a light-cone gauge propagator, 

r = f JlLpv(p±] l - (A 6) 

3 ' n -J(2ir) d ^(P+fe)((J-*)»+fe)((f-p)* + fe)' K '> 



where p + = pn, l + — n • I, and 'PV denotes use of the principal value prescription fl25 
in order to regularize the singularities of the light-cone gauge propagator. More precisely, 
since J^ n will be contracted with —g lJLl ,k v afterwards, we only need its part ~ fc M , as p, n do 
not possess non- vanishing components in their (d — 4)-dimensional parts. Assuming Lorentz 
covariance of the integral one finds 



■Kn = ~ ( Js,n -J3 + —Jl-n ) ^ + ■ ■ ■ , (A.7) 



2(1 - x) { J3 ' n J3+ \k 2 \ 

where J2. n , Js,n are the scalar two- and three-point functions with a light-cone gauge prop- 
agator as given in Eqs. (A.4) and (A. 14) of M, respectively, and J3 is the ordinary scalar 



three-point function of (A. 13) in g. The dots in (A.7) indicate the contributions propor- 
tional to p M and n M , which we do not need. We finally note that in our previous paper [Tl| ] 
on the polarized NLO splitting functions we were able to determine the results for most of 
the virtual graphs just from their contributions in the unpolarized case. The results in Ap- 
pendix A of m and the ones presented here now enable a straightforward direct calculation 
of all virtual contributions also for the polarized case. 



B Real Integrals 



The real emission processes p — * k + p% + P4 (see Fig. ||(c)) comprise integration over the 
momenta of two unobserved real massless particles. In addition to (A.l) we define 



P3 



(Pi 



Ps 



Pi , Ps ) 



(B. 
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introducing the (d — 4)-dimensional components of ps which we split into a part p 3 parallel 
to those of k and a transverse part p 3 . The phase space then reads (again including the 
integration over the 'observed' parton k): 



p 5 (2) 



d d kxS ( x 



pn 



dd P3 x +f 2^ 

^TT ( J (Pa) 



(2tt 



d d P4 



2 4 + 2e 7rr(l - 2e) 
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(27T) 



^ T «J+(^)(2 7 r) d 5 d (p-fc-p 3 -p4) 



x dk (k(l — «))" 



dkk 



dw (w(l — w)) 



dv (v(l — v)) 2 



where we have defined 



-3/2-e 



- / dA" (A^l-A 11 



(B.9) 



-1/2 



k 2 = \k 2 \(l - x)kk , 
k 2 = \k 2 \{l~x)k , 



P°3+Pl 



2P(1 - jc)iu = 2P 
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pn 
1 



{Plf-(PI? = 4+v(cl-cl) = -(p° 3 +pl)(p-p 3 ) 



p 



Pi 

±\2 



(Ps) 



Ai + A II (A 2 - Ai) , 
v(l-v) fa + afX 1 - 



with 



Cl,2 



Al.2 — 



\k 2 \(l-x)w 



XWK 



(B.10) 



((Pi) 2 ~ (Plf ~ c lC2 ) T ( Cl + c 2 )J(l-k)(l-X^)v(l-v) .(B.ll) 



Again the last three integrals in Eq. (BJ)) have been written in such a way that they are all 
unity if there is no dependence on (d— 4)-dimensional scalar products. If present, such terms 
only give contributions proportional to e after the last three integrals have been performed. 
The advantage of writing PS^ 2 ) in this way and performing the three last integrals first is 
that the result of these integrations can always be expressed in terms of usual (rf-dimensional) 
scalar products. For instance, 



1 

P 



p 2 ^ (Pi) 2 + (piy 



9fiu& Pz — k>P3 



1 - e 
2e 

1 - e 

e 



(\k 2 \ + 2x(p-k)) , 

(p-Pa) . 

pn 



(B.12) 



-(1 - x)(p-p 3 ) - 4- (p- k) + (p- k)— — 

1 — e V 2 pn 
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after integration over k, and A". The other two (d — 4)-dimensional scalar products 
that appear in the calculation, k 2 j)\ and (/cpj) 2 , have lengthy expressions which can be 



strai ghtforwardly obtained from (B.9). Thus, after performing the last three integrations in 



( B.9 ), there are only terms left in the matrix element which are familiar from the unpolarized 



case. The matrix element can then be further integrated using the results of Appendix B of 



or, equivalently, doing the remaining integrals in Eq. ( B.9 ) 
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